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THE PARTIALLY AVERAGED FIELD APPROACH 
TO COSMIC RAY DIFFUSION 

Abstract 

The kinetic equation for particles interacting with turbulent fluctuations is 
derived by a new non-linear technique which successfully corrects the difficulties 
associated with quasi-linear theory. In this new method the effects of the fluc- 
tuations are evaluated along particle orbits which themselves include the effects 
of a statistically averaged subset of the possible configurations of the turbulence. 
The new method is illustrated by calculating the pitch angle diffusion coefficient 
D^^ for particles interacting with "slab model" magnetic turbulence, i.e., magnetic 
fluctuations linearly polarized transverse to a mean magnetic field <B> . Results 
are compared with those of quasi-linear theory and also with those of Monte Carlo 
calculations reported in a companion paper. The major effect of the non-linear 
treatment in this illustration is the determination of D^^ in the vicinity of 90® 
pitch angles where quasi-linear theory breaks down. The spatial diffusion coef- 
ficient K|| parallel to <B> is evaluated using D^^^ as calculated by our technique. 
It is argued that the partially averaged field method is not limited to small ampli- 
tude fluctuating fields and is hence not a perturbation theory. 
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THE PARTIALLY AVERAGED FIELD APPROACH 
TO COSMIC RAY DIFFUSION 


I. Introduction 

The systematic treatment of the transport of charged particles in disordered 
magnetic fields began about ten years ago.‘*^ These authors derived a Fokker- 
Planck equation for the time evolution of the particle distribution function by a 
method that is commonly referred toas the quasi-linear method.'* 

A more formal definition of the quasi-linear approach will be given in the next 
section but the basic idea may be quickly sketched. To compute the mean square 
effect of a random force field on a particle one must be able to follow the particle's 
orbit in that field for such a time as is required for the random force to become 
self-incoherent. This is an inherently non-linear problem and in general intract- 
able; however, if the random force can be separated into an average part, for which 
orbits can be calculated, and a random part that is small (in some sense) one may 
use the orbit in the average field to compute the effects of the random part in the 
hope that the random part will become self-incohereni before it has a chance to 
significantly perturb the orbit. 

There are two approaches one may take in applying these notions. One is to 
assume that the correct equation of evolution is a Fokker-Planck equation and 
proceed to calculate the appropriate coefficients. This approach has been de- 
veloped primarily by Jokipii*'® and by Hasselmann and Wibberenz.^ The diffi- 
culty in such an approach is the proper treatment of non- random forces.^ The 
second method is to begin with the Vlasov equation and derive a Fokker-Planck 
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equation to describe the evolution of the ensemble averaged distribution function. 
The ensemble is the collection of all possible random fields with the required 
statistical parameters. This approach was employed by Hall and Sturrock, ^ 
Roelof,^ Kulsrud and Pearce,® and discussed in a general analysis by Kaufman.® 
Most recent authors have employed this approach and we shall use it in this paper. 

The first questions about the appropriateness of a Fokker-Planck equation for 
describing charged particle transport in random magnetic fields were raised by 
Klimas and Sandri.*'' They pointed out that the presence of particle trajectories 
(in the average field) with pitch angles near or at 90° with respect to the average 
field invalidated the "adiabatic hypothesis." Since this hypothesis is necessary 
in order to derive the Fokker-Planck equation from a non-Markovian integro- 
differential equation resulting from previous steps in the derivation, these authors 
asserted that the correct transport equation was the more primitive (and more 
complicated) non-Markovian one. 

The present authors then pointed out” that these same 90° pitch angle particles 
caused even worse problems. Such particles remain in a region of correlated field 
for an arbitrarily long time and thus their orbit perturbations could become im- 
portant before the random force becomes self-incoherent, violating the assump- 
tionsofquasi-lincar theory. Indeed we were able to show that for a particular, 
simple model the terms that are dropped inquasi-lincar theory eventually domi- 
nate for such particles. 

Although the difficulty with quasi -linear theory is confined to a region of phase 
space that is usually small, this region about 90° pitch angles can be quite im- 
portant in considering the bulk transport of the charged particles. Earl’^ has 
shown that if the fluctuations in the magnetic field are depleted at high frequencies 
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(power spectrum steeper than k~^) quasi -linear theory predicts no scattering 
through 90° pitch angles and hence a partially coherent propagation of the par- 
ticles rather than a diffusive propagation. 

Thatquasi-linear theory is inadequate to describe the behavior of particles 
with pitch angles near 90° is generally accepted and several authors have pro- 
posed modifications designed to obviate this problem. Some of the proposals 
have been purely phenomenological,*^ some formal*'* and some have been mix- 
tures of the two.*^ It is the purpose of the present paper to describe and pre- 
sent the results of a new method of calculating the diffusion coefficients of 
charged particles in a random magnetic field. 

This method, which we call the partially averaged field (PAF) theory, has 
been briefly described in a previous report.*® Its fundamental idea may be 
summarized quite simply: while in quasi -linear theory the effects of the fluc- 
tuating field are computed along a particle orbit in the ensemble averaged field, 
in the present approach the orbit that is used is the particle's orbit in a partially 
ensemble averaged field. The partially averaged field is averaged only over a 
subset of the full ensemble, namely the subset of all realizations in which the 
value of the field at a certain fixed point has a fixed value. The results com- 
puted using this orbit are then finally averaged over the rest of the ensemble, 
i.e., all possible values of the field at the fixed point. 

The resulting theory is non-linear in that effects of the fluctuating field are 
included in the particle's orbit in a fundamentally non-linear way and has the 
virtue of satisfying, at least marginally, the requirements for a useful, approx- 
imate theory. 
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These remarks will be made more precise in the next section where we will 
discuss the formalism employed in deriving the kinetic equations and show how 
the partially averaged field theory fits into this formalism. The conditions under 
which phase space diffusion takes place are developed. As an example of the 
application of our formalism we consider in Section III the "slab" model of the 
fluctuating field. lATiere existent, the pitch angle diffusion coefficient is explicitly 
evaluated and compared with the result of quasi-linear theory. The differences 
between the results of quasi-linear theory and partially averaged field theory are 
discussed in terms of the physical processes involved. In Section IV we review 
the relationship between pitch angle diffusion and spatial transport. For the pitch 
angle diffusion coefficients derived via the PAF theory in Section III, we derive 
in this section the corresponding spatial diffusion coefficients parallel to the 
mean magnetic field. Since depends importantly on the pitch angle diffusion 
coefficient in the vicinity of 90° pitch angles, W'e develop in this region the approx- 
imate scaling (from PAF theory) of the pitch angle diffusion coefficient with both 
the rms level of magnetic field fluctuations and with particle rigidity. Finally in 
Section V we argue that PAF theory is unrestricted in its validity by the rms 
value of the magnetic fluctuations and is — at least marginally — a genuine non- 
linear theory. 

II. Derivation of the Kinetic Equation 

Consider a collection of charged particles of rest mass m, charge q, and 
momentum p [/ » (1 + pVm^c^)*-^^] which move in response to the electromag- 
netic force F * q(E + v x B//mc). In our later cosmic ray considerations the 
electric force is negligible compared to the magnetic; for the time being, however. 


4 







we proceed with generality. F has a component which fluctuates randomly in 
space and/or time. For our cosmic ray application of the theory F is thought 
to be an "externally imposed" force, although the additional constraint of Maxwell’s 
Equations can be added straightforwardly. 

We describe the particles by their one partic’e distribution function f(r, p, t) 
which, when the particles are interacting only with an externally imposed F, obeys 
exactly the Vlasov Equation 



( 1 ) 


We are interested in the coarse-grained behavior of f - its evolution on the 
time scale of several coherence times between the particles and the fluctuating 
force. We invoke the usual ergodic hypothesis and assert that this behavior is 
identically the behavior of ( f > , which is f averaged over an ensemble of F's. 
The ensemble of F’s is specified by the values of the moments 


<F) 


<FF) ^ '((F) 4 .-f) (<F) +bF) - <F> <F> + <bFoF> 

(FFF) (F) <F> <F> + (F) <bFbF> + (bF <F) bF) -t C-FoF) <F) 

j r\j ''j ^ ^ ^ 

-t <jFcF:;F) 


etc. 


Equation (1) is readily split into an equation for : f ) 


{— -V + <F> • — <f> 

V 3 1 my ~ bp y 



b f 





< a F 0 {) 

( 2 ) 
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( d/,^ p . F * 0 for the electromagnetic force) and one for cf »= f - ( f > 



<F) 





Equations (2) and (3) together contain the full information content of (1). 

An equation for < f > alone is obtained if Eq. (3) can be solved for f in terms 

of Its initial value f(r, p,t^ )and <f) and the result substituted on the right side 
of Eq. (2). An exact solution for f is impossible, for Eq. (3) mixes derivative and 
ensemble average operations in a non-commutable fashion. 


(a) Quasi-linear theory 

The quasi-linear method approximates S f by the solution to the equation 


(F) 




( 4 ) 


which neglects terms in Eq. (3) bi-linear in the fluctuations. The corresponding 

equation for /f ) is 

\ / 




The quasi-linear propagator U^lt.t') propagates the r. p dependence of functions 
to which it is applied backward in time from t to f along trajectories in the force 
field ^F) ; these trajectories are the characteristic trajectories of Eq. (4). 





We discard the initial value term from further consideration on the assump- 
tion that after a reasonable interval t-t^ all correlation between sF and U^lt, t^,) 
if(r, p, tg) has died away so that ensemble averaging yields 0. If such correla- 
tions do in fact persist, we think that it is impossible to describe < f > by the 
simple Markovian description toward which we strive. We assume that a similar 
argument applies when Ug is replaced by the more correct propagator that we 
introduce later. 

It is the spirit of quasi-linear theory that scattering is small over a partic le- 
fluctuation coherence time . The evolution of < f > in Eq. (5) is predominantly a 
convection of ( f ) (tg) along the characteristic trajectories in < F) . The quasi- 
linear collision term 


COL 


-• I <1 • U (t, r)cF(r, p, T) • — <f> (r. p, 

P ' - . 'p - 


. 


( 6 ) 


represents a small correction. It follows then that the adiabatic approximation 
Ug(t, 0 Mf. P. - < f ) (r, p, t) can and should be made in every instance that 
the quasi-linear theory is valid.*** The quasi-linear collision term is thus prop- 
erly written 


C*?*- r _ 


'I'- (t. T) P* • — ><f) fr. p, t)(7) 


P 


The brackets i ^ indicate that Ug now' propagates F • / 'p only. 

We now examine the circumstances under which quasi-linear theory is valid, 
i.e. is "small." It can be shown'* that quasi-linear theory is valid for most 
regions of the particle phase space provided that E and ,B are small in the 
sense ''efined in Ref. 18). The reason is that the correlation \ F(r, p, t) U^(t,') 
P*^)/ diminishes rapidly with t - r. This decrease occurs in the 
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characteristic time t^' and may result from (1) a short intrinsic correlation 
time, i.e., the correlation < F (r, p, t; '*F (r, p, t)) at fixed r, p cuts off rapidly 
with increasing t- t, and/or (2) a short intrinsic correlation distance in r through 
which the quasi-linear characteristic orbits move rapidly in time. The condition 
for the validity of quasi-linear theory may be represented concisely as t^*- <■' 
t where t is the characteristic scattering time for particles in the pres- 

sea t ’ scat 

ence of the fluctuating force. 'WTien t«L « t^^^^ , scattering effects can be 
neglected in the computation of the scattering itself; quasi-linear theory is an 
iterative procedure. 

The quasi-linear collision term Eq. (7) can be further simplified by noting 

that 


where p( ? ) 


Uo(t, Op *6 the propagated momentum. 


’pr O 


r(T ),1 


P . J1 ^ 

'P(t) .^p 'p(r) 


Further 
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since a variation in p( ' ) at fixed r( ) results in a change in both r and p. The 


quasi-linear collision term is in general thus composed of two components 

r 

01 si 1 > 

(•FUod. J) p. r)^ 

. X . 

< f ) ( r . p, t ) 

CvJ- - • 1 (J f < 

t 1 

..lp(T) 

P 

r' 

k It 

<^'F Eq (t. r) 1 F (r. p. t)''^ 

'r 

s • 

- ( f > (r. p, t ) 

4 • 1 UT 

T- J.o 

^'p(0 

"r 
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The first term represents pure momentum diffuc’on. For times t-t^ » t^*- 
this diffusion is described by the tensor 


DQl 

pp 




p, t ) U ( t , t - 7 ) { ;- F ( r , p. t - T ) 


'> 




(9a) 


There is in addition a mixed p, v diffusion, which for large t-t^ is described by 
the tensor 


DQL - 

r.pr 



•F (r, p, t ) Uq (t, t - 7 ) • ■ F (r, p, 



(9b) 


This mixed p , r diffusion disappears when < f > is spatially homogeneous and 
may be small for weak inhornogeneities. It is generally present, however, to the 
same degree as the pure momentum diffusion. 


(b) Partiall y averaged f ield theory 

In many models there are, however, pathological though important regions of 
phase space where rj:’* is quite long, so long in fact that ‘j” t and 

both the spirit and the substance of quasi-linear theory are ’ iolated. 

We propose here an improved approximation, which includes the essential 
effect of the fluctuations themselves in the lowest order propagator. For reasons 
that will be soon evident, we call our method the partially averaged field approxi- 
mation. In this formalism the coherence time j is naturally limited ^ 
and the collision integral is validly calculated for (1) all regions of the particle 
phase space and (2) arbitrary amplitude of the fluctuations. Our method is not a 
weak-coupling approximation. It reduces to quasi-linear theory when r[:^‘ ^ 

holds. 
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Our formalism proceeds as follows. We consider the ensemble averaging 
process to occur in two steps. A first partial average, denoted by ( >y , is per- 
formed over all realizations with the same fixed values of S E(r,t) and h B(r,t) at 
the point r and time t at which Eq. (2) is valid, A second final average, indi- 
dated by < ) ; , is then carried out over the distributions of ^E(r,t) and 6B(r,t). 
The collision term in Eq. (2) is identically written according to this prescription 


C 


^ • (bF(r, p. t) ( ~A) n 

dp 


( 10 ) 


If if(r,p,t) depended only on dF(r,p,t) it would be the same for every real- 
ization included in the partial ensemble and partial averaging would be trivial. 
However, f depends on oF globally. We therefore write Eq. (3) for the arbi- 
trary point r', p’, t' and partially average to obtain 

\ 

7^ ' A 

< i) f > ._ + >>F ’ of V - ' 'F • fc f > 


P n' /E-N 

+ 'V + ( F > 


-t' 


nr>' 


■*p 


= - < T'F) 



( 11 ) 


Note that partially averaged quantities depend functionally on r', p’, and t' which 
are the variables of Eq. (11) and parametrically on r and t which are singled out 
by the partial averaging procedure. Thus far Eqs. (10) and (11) are exact, for all 
we have done is to break the full ensemble average procedure into tw'o sequential 
steps: symbolically < ) => <^< > ' i 

The essence of our method is to approximate in Eq. (11) as follows 


/ 7. '\ i > \ 

( oF • bf) , - / <■ F • — f f ) ^ (iF) . • — (if) 

\ r.p' / \ '■' ip' / ^ Tip' 


( 12 ) 
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This approximation maintains the properties of the exact terms in two essential 
limits. First, when r*. t' is much different from r, t partial averaging and total 
averaging are equivalent and < 3 F • 9/7)p* 3 f > - < 3 F • ^/^p' 3 f > * 0; however , 

< F\, = 0 also under such circumstances. Second, at r' * r, t' * t we expect 
<S F • S/r)p' 3f ) j (- S F • 3/9p* <3 f ) to dominate < 3 F • rf/Bp* 3 f > , particularly 
for those realizations that contribute importantly in Eq. (10); the right side of Eq. 
(12) is, however, identically 3F • B/^p' <*f > ,,, at r' *= r, t' * t. 

Our approximate equation for <S f > ,,(r', p', t*; r, t) is thus 

I ^ 



dt' my' 


V' + (<F) 

' 'V 


\ 

+ (5F)r,j • — 
' Bp' 


B (f) 

<3f) p = - <SF) p • 

^ Bp' 


(13) 


Neglecting the initial value term, whose contribution to C is again assumed to 
damp rapidly with time, we express the solution to Eq. (13) as 


< j f > .^ (r', p', t' ; r, t) 


f 


dTU.;(t',T, t) < 3F> p fr', p', t; r, t) 


B (f> 

(r', p', r) (14) 

Bp' 

The operator U. (f, t; t) is the propagator in the partially averaged field. 
Operating on an arbitrary function g(r', p', t) it propagates the r', p' dependence 
backward in time from t' to r along orbits in the partially averaged force 
< F > + < F ) . U also '’epends parametrically on t because the partial average 
is carried out at constant (■ E(r,t) and 3 B(r,t). 

For Eq. (10) we need < f (r' » r, p* * p, t' - t; r,t). Substituting from (13), 
our new collision term becomes 
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We expect decay of <(s F U^j (t, 7 ; t) < S F > p (r , p ,t ; r , t)^ in the time t - t » t^* 
(tp* 1 t^^) sufficiently short that all characteristic orbits in the partially aver- 
aged force remain at least marginally valid for the entire time interval over 
which the integrand in (15) is non-zero. We hold this expectation irrespective of 
the size of SF and shall argue why in Section V. 


For times t - t ^ we again write 



The adiabatic approximation on < f > cannot always be made in C . The physi- 
cal reason is that when 5 F is substantial, scattering as described by Eq. (16) repre- 
sents a significant correction to the convective motion in < F ) . Our success in de- 
riving an equation unrestricted by the size of SF is thus somewhat mitigated by 
the fact that this equation for < f > is an integro-differential equation. One must 
not conclude, however, that an adiabatic approximation is never possible: each 
model to which the < f > equation is applied must be examined individually to see 
whether an adiabatic ansatz can be made consistently with the resulting physics. 
Work reported in later sections is illustrative of this point. 
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(c) Gaussian statistics 




It is impossible to proceed further with the partially averaged field method 
in full generality. We always consider the fluctuating electric and magnetic fields, 
& E(r,t) and h B(r,t), to be Gaussian processes. It can be shown (Appendix A) for 
a Gaussian process that 


< SE ) (r', t'; r, t) = &E (r, t) • Cj. (jr', 


(r', t'; r, t) = iiB(r, t) • C„(r', 


t', r. t) 

(17a) 

t', r, t) 

(17b) 


where C,, and C„ are the normalized correlation tensors for the fluctuating 
electric and magnetic fields respectively. 

The assumption of a Gaussian process is a common one, partly because the 
central limit theorem*^ states that a tendency to\\'ards Gaussian statistics is a 
common occurrence in random systems, but mostly because of its mathematical 
tractability which most other types of statistics do not share. Mathematically, 
Eqs. (17) follow uniquely - and bE and bB are Gaussian processes - provided 
that one demands that < S E > (r*, t': r, t) and < b B > (r',f; r,t) be linearly de- 

pendent respectively on feE(r.t) and b B(r,t), 


III. The Magnetic Slab Model 

We have worked out the details of the partially averaged field formalism for 
the so-called magnetic slab model, and we shall discuss details in this section. 
The magnetic slab model was chosen because its simplicity facilitates calculation 
and direct comparison with Monte Carlo results (which are reported in the ac- 
companying paper and because it also contains some of the basic physics of 
the scattering of cosmic rays from irregularities in the interplanetary medium. 
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Tn the magnetic slab model there are no electric fields. The ensemble aver- 
age magnetic field < B > is spatially homogeneous and oriented along the z-axis of 
a Cartesian coordinate system; magnetic fluctuations are linearly polarized, con- 
stant in time, and vary in magnitude only with distance along <B) . Thus 
SB* oB (z) e^. In this simplest form, there is no spatial variation in planar slabs 
oriented perpendicular to the z-axis and hence the coinage "slab model." At the 
end of this section we discuss a more generalized slab model in which S B is in 
the e'^ direction but planes of constant | SBj are! to the y-z plane and at a fixed 
angle (tjpically 30°) with respect to the z-axis. 

The magnetic slab model is intended as a description of the interaction of 
cosmic rays with slow magnetic structures (tangential discontinuities or Alfven 
waves) in the solar wind. Let us assume that the radial divergence of the solar 
wind flow can be neglected and that we can place ourselves therefore in a unique 
solar wind frame of reference. The amplitude of the electric fields associated 
with these slow structures is at most i E | = v^ 1 B | /c and hence far smaller than 
the magnitude of the magnetic fields (v^ is the Alfven speed and is tjpically 50- 
100 km/sec). The magnetic force on relativistic cosmic rays is therefore far 
larger than the electric force, and for that reason we neglect the latter. 

Our assumption that <B> is completely uniform breaks down, of course, on 
the 1 AU spatial scale of the classical Archimedes spiral. However, the pertinent 
scale for a description of scattering is the autocorrelation distance of magnetic 
fluctuations, which has been measured to be a few hundredths of an AU;^* to the 
extent that < B ) is uniform ever this length our model is a good one. There is 
also evidence that the magnetic fluctuations are preferentially polarized in the 
direction normal to the ecliptic plane (the normal direction would thus correspond 
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to our x-direction); any component of bB in the ecliptic plane and normal to <B> 
at its point of generation gradually develops a compressional component as it 
propagates outward in the solar wind and tends to damp away.^^ 

Since electric fields are negligible, cosmic rays suffer no energy changes in 
their interactions with the fluctuations and |pl is simply a parameter in the equa- 
tion for <f ) . Scattering in angle does occur and can be described in terms of the 
variables P and <p= tan~‘p /p . is the cosine of the cosmic ray's 

pitch angle relative to < B> and t is the gyrophase angle measured from cB (the 
x-axis). 

We consider here the situation where (f ) is spatially homogeneous. Such a 
situation arises when ( f > is homogeneous at t «= t^ and the spectrum of the fluc- 
tuations is itself spatially homogeneous. This assumption of spatial homogeneity 
is artificial and unrealistic: as a result of it the bulk spatial transport of cosmic 
rays along and across field lines is precluded. It is made solely to achieve a 
simple model in which to explore and justify the partially averaged field method. 
The spatially inhomogeneous situation will be discussed in Section IV. 

With the aforementioned approximations the equation for (f> is 

Mf> ' (f) /’ r’ 

- < . > . • ■ D (/., /) < •, V (It U (t, t - r) b , 

\ Ji) 

\ 

Dir, f) <f> (r, p, t-T)'» . (18) 

■■ - / 

In Eq. (18), ( . > q(B) /, me is the gyrofrequency in the ensemble average mag- 
netic field <B)e^, is the gyrofrequency q j. By/mc that would prevail if 
I Bie^ were the only magnetic field, and D is the differential operator 
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D = (1 _ MCOS ^ 


() 

dfX 


(1 sin 


9 fi cos t 

(1 


(19) 


We next approximate further by neglecting the gyrophase dependence of < f ) . 
This is accomplished formally by first expanding < f > in the Fourier series 


<f > 



< f > ^ exp in//, 


next substituting into Eq, (18) and projecting out the n = 0 component 


^ < f> 

(3t 


r 2 t 7 


dfi 


(1 - sin 


X.,. 


•'0 


dr U;,^ (t, t - t; t) 



2\l/2 


sin 'f. 


r) 

3/i 


COS 

(T~~^ 




(f> 


. \ 

exp 1 n /') . 

n ' ' j ' 

/ 


(20) 


and finally neglecting the n 0 terms on the right hand side of Eq. (21). Fourier 
amplitudes < f > ^ , n > 0, are small compared to < f ) q in the approximate ratio 
■= < i B ^ >/( B) ^ (cf. Appendix B) provided that they are initially small; the 
physical reason being that the convective term < ^ > :>< f keeps particles well 
stirred in / . Thus 
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dt 


J_ 

2t7 



— (1 sin /. 

■-/X 


/ r 

( du,> •J''" Upj (t, t - T ; t ) 

\ ^0 


c ' (1 -/X sin / • 


<f> 

dii 


0- ^)>:i 


( 21 ) 


Next note 


so that 


-4 


(f ) 


dt 



( 22 ) 


(23) 


To the right de of Eq. (23) must be added any effects such as sinks and 
sources of particles which are not included in the original Vlasov description 
Eq. (1). 

In the steady state < f ) ^ no longer depends explicitly on time and (23) can 
therefore be written in the more transparent form 


Here 



f 2rr 


- 

1 

277 

d/. <( 

: <f> .) ) 



0 



\ 


’ <f' oOO 



(24) 


d ( f > 


dr 


dr 


is the change in < f ) g apparent to an observer moving from the point p at t * 0 
along an orbit in the partially averaged field < ?B) . 

To the extent that < f > ^ (/j) is a linear function of ,, - and we have found from 
our Monte Carlo calculations that such is approximately the case for a wide variety 
of conditions - 
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0 


<^A (f > <A/;>,, 


d/2 


and Eq. (24) is of diffusion form 

r> <f> 


d/i 


DPA 




1 

— A 

1 




2rr 

(a <A;,).^ 


d ( f ) 


( 25 ) 


The diffusion coefficient D has exactly the same form as the quasi-linear 
diffusion coefficient but differs in that A/j is the change in evaluated fol- 
lowing an orbit inthe partially averaged field rather than following a helical orbit 
in (B) as quasi-linear theory prescribes. For a /^ = 0 particle in the slab model 
A;. approaches no finite limit as t - oc in quasi-linear theory^^ because such a 
particle suffers exactly canceling deflections over its circular orbit at constant 
z in ( B ) ; no such pathological behavior occurs in the partially averaged field 
method because even /i ■ 0 particles move freely in z and hence enter regions 
of uncorrelated 5 B. 


(a) Results for the Perpendicular Slab Model 

In order to evaluate < A,; ) . and hence the diffusion coefficient 


DPA 
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Newton's Equations of motion 


^ - p' X [< . > ^ . (z) C (z - z' ) ] 


(27a) 


<lr’ p' 
(It /tt) 


(27b) 
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are integrated numerically backward in time from the boundary conditions 
P'(t ss t) » p, r'(T » t) » r. The correlation function C is chosen to be the 

<*w 

exponential 


C (z - z' ) 


1 

exp - 

z 



(28) 


since we compare with the results of a Monte Carlo simulation,^" for which 
it i." simplest to generate numerically magnetic fields with an exponential corre- 
lation function. The correlation distance z is the distance scale and the corre- 

C 

lation time t * z > m/p the unit of time. 

We first pick a b l (z) and the boundary values * t) and / ( r * t). / for 
inclusion in (26) follows immediately from Eq. (22). < ) , on the other hand is 

the change in /i which results from the integration of Eqs. (27). This integration 
is carried out until the particle has moved tjpically a distance z' - z ^ 16z^. in 
the direction of varying o - (z’). 

The entire procedure is next averaged over a uniform distribution of ; (t). 

The density w'ith which the ; (t) must be chosen depends on , (t). WTien the latter 
is small (/. (t) X (z)/< - )) particle mirroring is an important effect, as we shall 
see shortly. For small ,, , the ;-*s are taken at 1° intervals in order to achieve 
convergoiice for the /-averaging process. For larger (X (z)/< < )) mirroring 
does not occur and it is sufficient to take /(t)'s w'hich are 2° apart. 

The last step in the evaluation of Dj”* is the final averaging procedure, in 
which we perform the steps of the preceding two paragraphs for several different 
values of . (z), typically 20, each weighted by the distrilxition 

1 • . ^ 

P( .)■ <‘xp 
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Figures 1 to 3 depict thus computed as a function of Each figure 
corresponds to particular values of 7 / * < > and ► ■ > m ( - > z^p. 

The discrete points at which is calculated are indicated. The solid curve 
is a spline fit to these points. Shown also on each graph is the curve of D vs 
fj. This quasi-linear diffusion coefficient is computed from the formula 


D<?l 


1 -/ 


f. ^ -t 


(29) 


which results when the exponential correlation function (28) is assumed. 
Equation (29) is the limiting t • ■ value of the expression 


wherever that limit exists. As we have mentioned earlier and demonstrated else- 
where this limit does not exist as/. - 0 and indeed a quasi-linear diffusion co- 
efficient is not properly defined at such small /i’s. 

We would point out three aspects of these figures: (1) the peaking and definitely 
non-zero value of as ,, - 0 . ( 2 ) the fact that is approximately equal to 

but slightly less than D^,'; for • , and (3) the essential equality of D'*'' and 
as /■ 1. These aspects are also characteristic of the diffusion coefficient 

obtained from Monte Carlo computer experiments. 


(b) Interpretation of the Differences between PAF and Quasi-Linear D 's. 

— un 

We interpret the bump which appears in the vicinity of , » 0 as due to par- 
ticle mirroring: particles with pitch angles near 90® are dominantly affected by 
the maximum in ( B> at z * 0 and their / 's are determined by magnetic moment 


conservation to be 


I 


I 


i 


< A/^ ) ^ 1 

(Without loss of generality we assume here that z = 0 is the point at which the 
diffusion coefficient is to be evaluated.) Here »= oB(0)/ (O and is the co- 
sine of the particle's pitch angle with respt 't to the total (partially averaged) 
field ( B > 
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» 2 


(UV) 


1/2 




- M 


(30) 
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cos A 


( 31 ) 


Using Eq. (30) for A/x and Eq. (22) for /' , we determine the mirroring contribution 
to D|^ thus to be (cf. Eq. (26)) 



The T sign is included within the integral sign because the algebraic sign of <p> ■ 
is dependent on /; if it were not, (32) w'ould clearly integrate to zero because n' 
depends only on cos / , and there would be no contribution of mirroring to the 
diffusion coefficient. Furthermore the sign cannot be just a function of cos : (as is 
the sign of /x' ) or the same reasoning gives a zero result once again. We shall see 
that it is just the act of mirroring that injects the needed asymmetry about * 0, 
i.e., a dependence on sin to produce the peak of tJ » 0. 

To determine the sign of < /x ) . we proceed in the following heuristic fashion. 
Consider the /. > 0 (but | ^ 0) particle at z - 0 at t * 0 in the partially averaged 
field. Depending on its ; and /x, this particle's guiding center is at positive or 
negative z 
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(33) 


and its ,/ (cf. Eq. (31)> is also positive or negative. If i. „ and are both > 0, 
its initial motion carries the particle to the right (to larger z) and we assert that 
</. ) , > 0. Similarly, if and ,/ are both ) •' 0. On the other hand, if 

z and // have opposite algebraic signs, the particle's initial motion carries 
its guiding center toward z = 0 and stronger magnetic field so that mirroring is 
a possibility. If mirroring does occur — and we shall presently develop criteria 
for mirroring -(,.)■ has the same sign as z^^. If mirroring does not occur, 

< /. ) . has the same sign as . In Figure 4, which depicts the transverse p__ , p^, 
plane, mirroring may occur for ) > 0 in those sectors denoted by an "M?". 
(Because of the symmetry of the distribution of B(0), we may restrict our con- 
sideration to , > 0.) The sign of ' /. > in the complementary sectors is also 
indicated. The angle A is given by 


Sin 




2x1/2 


When does mirroring in fact occur? We first applied the criterion that pre- 
dicts mirroring when 


1 > 1 

<B> < B) . (z 0) ' 

On using the approximation 

(B) (z^-c ) - <B> (z- 0) 



(34) 
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plus the expressions (iJl) nnd (33) for /j' and the criterion for mirroring 

is equivalently written 




2^l/2 


2 ^ 1/2 


COS 


n -/d) 
(1 ^V) 


\ Sh' 


(35) 


The : -integration in Eq. (32) has been performed numerically. The ambiguity in 
the sign of <,, > in the regions of Figure 4 r rked by "M?" is resolved by test- 
ing for the validity of (35) and assigning (. ) the sign of z if (35) is indeed 
satisfied. 

Following the procedure just prescribed, we have calculated the mirroring 
contribution to Dj;;* and found it in the vicinity of « 0 to be typically a factor of 
2 larger than the determined by exact integration of orbits in the partially 
averaged magnetic field: mirroring has been introduced too strongly in our 
heuristic procedure. Much better agreement is obtained when we modify the mir- 
roring requirement (34) as follows: 


' > \ (30) 

The additional factor in the denominator of the right side accounts for the fact 
that even when the guiding center is located at z * 0 the particle itself is almost 
ahvays at a weaker field strength, owing to the exp - ( /. /z^. ) fall-off in D . 

We attribute the fact that D0‘- in the regions > of Figs. 1-3 to 

another effect, which we shall try to explain heuristically in this section. In 
quasi-linear theory, a particle moves with the assumed constant momentum p,. 
in the z-direction. As time increases from t =* 0, the particle interacts (in the 






<B> (z 0) 
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ensemble average) with a portion of the power spectrum of the fluctuations of 
ever-decreasing width: the effect on of all Fourier components except that 
at the resonant wave number k = e/i/^ | z ^ ultimately averages to zero. 

The constancy of /i and p^ are, however, only approximate. More accurately, 
the particle attempts to maintain the constancy of /^' , Eq. (31), which we write in 
a slightly modified form as 

— ~ ^ (1 ’) cos /.' = const. (37i 

p d T 


Equation (37) can be integrated approximately to give 


z = /.; + (1-/; 2)1/2 cos .'I T _ ^ ~ ^ '' [sin(/ + </> T) 


- sin •/.] (38) 

Note that this z-motion includes an oscillation at the gyrofrequency < - ) , As a 
result, the particle never attains perfect resonance with a single Fourier mode 
and its interaction with the fluctuations is on the average diminished. 

To see this semi-quantitatively, we note that by approximately integrating 
Eq. (22) A/j is given as 
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Equation (39) differs from its quasi-linear analogue only in the fact that the oscil- 
latory component of the z-motion is included, ';/in (38) has been approximated by 
t] and the -q contribution to the rectilinear motion has been discarded. 

It is now straightforward to calculate, using Eq. (26), the resonance contribu- 
tion to We use (22) for /' , (39) for ' /^, assume no correlation between dif- 

ferent Fourier modes (the random phase approximation), and employ the expansion 


exp (- i 


+OT 



m Cl 


exp - i m ,■> (a) 


where the J^^'s are the ordinary Bessel Functions, The result obtained for the 
exponential correlation function Eq. (28) is 


pr.-s 

Ul-i 
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(40) 


The quasi-linear result is recovered by letting q -• 0 inside of the summation, only 
the s * 1 term survives and Eq. (29) readily follows. 

It is not obvious from Eq. (40) that the resonance interruption results in a 
diminished diffusion coefficient. We have, however, evaluated numerically as a 
function of /. the ratio 
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for several values of ;; and e . Its value is indeed and its general behavior is 
quite similar to that of in Figs. 1-3 in the ,■ > // range. This finding is 

MM M M 

PA 

the basis of our interpretation that resonance interruption is the reason that 
falls beneath in this range. 

(c) Results for the Oblique Slab Model 

We have also evaluated D using Eqs. (22), (26), and (27) for the situation 
where 

r: bB(s) C(s -s') (41) 

The correlation C is again exponential (cf. Eq. (28)) but in this instance 

s - z cos a + y sin fx 

SO that planes of constant I (' b) 1, are perpendicular to the y-z plane and inclined 
at angle ' with respect to the z-axis. (The perpendicular slab model considered in 
Ilia corresponds to j = 0). Since the final averaging is over an assumed spatially 
homogeneous distribution of (s), the position s is again irrelevant and we may 
without loss of generality in (41) take z *= y » s = 0. The processes of averaging 
numerically over ; and oB(s) are exactly as detailed in Ilia. 

We display in Figure 5 results for and for d*e case t * 1,0, = 0.1, 

and r: » 30°. Figure 5 is thus the counterpart of Figure 2 for the oblique situation. 
The obvious difference betw'een the oblique and pc ’jjandicular models is the strong 
suppression of diffusion at small ^'s which non-perpendicularity introduces. This 
is true for both and . In additioi; -perpendicularity apparently leads 

to enhancements of in the region of maximum around f. » 0.5. The discrete 
values of /i for which D^^'s were explicitly evaluated are indicated in Fig. 5. The 
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distribution and number of such points is dictated by the length of the numerics 
involved in calculating For comparison we have also calculated at the 

same values. It is probable that the relatively small number of points in the 
region /j = 0 leads to the suppression of snr all amplitude oscillations which char- 
acterize the higher harmonic (of < ■ > ) interaction of particles and magnetic per- 
turbations with spatial structure perpendicular to . It is known that such 
oscillations characterize and it is probable that they would also appear 

if were calculated sufficiently densely in /x. Our attitude is that if these 
oscillations have the same amplitude in as they have in they are in- 
consequential by comparison with the magnitude of {/x ^ 0). 

As we shall see in Section IV, the spatial diffusion coefficient ^i, parallel to 
(B'^ varies inversely with . The important inference to be gained therefore 
from Fig. 5 is that the oblique slab situation is more conducive to parallel spatial 
transport than is the perpendicular slab. Extrapolating to the case of the inter- 
planetary medium, we expect that quasi-planar magnetic structures propagating 
radially outward from the sun undergo transition ‘^rom a perpendicular slab-like 
situation to the oblique slab case as the idealized <B> develops spirally. Given 
constant ■ and t should therefore ideally decrease with approaching distance 
to the sun. 

IV. Spatial Transport and the Scaling of (/- ^ 0) 

Let us assume that pitch angle diffusion theory is true locally and that 
is properly given by partially averaged fie.d theory. We re-introduce a spatial 
gradient in the z-direction with scale much longer than , and suppressing the 
suffix "0" rewrite Eq. (23) as 
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Contact with the macroscopic variables - density n, flux «J> etc. - is made by 


taking moments of Eq. (42) with respect to m. Particle conservation 
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(43) 


follows exactly by taking the zeroth moment and noting that there can be no dif- 
fusive flux 3 a > /'dfi at the boundaries /i * ±1 in m space. 

There exist a variety of methods for expressing <Ij = d^/ <f > in terms of 

2 

n and thus obtaining a closed equation for n. The standard transport theory method 
is to expand <f> in Legendre Polynomials and truncate the series after Pj (m). This 
yields a standard diffusion equation for n 


where 
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Alternatively one could expand <f > in eigenfunctions of the diffusion operator 
D rt/ Mu and obtain an alternate expression for m.? 
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where the t. are the eigenvalues of the equation 


and*^ 
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(48) 


Due to the difficulty in obtaining sol itions to (47), however, the usefulness of this 
approach is limited. 

A third approach, due to Jokipii,* is to solve the approximate equation 


— <f> = — D — <f> 

^ Bz B/^ 


(49) 


by inverting the diffusion operator to obtain 


<f> 


= f 

•'n 




(50) 


where we further approximate 
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Inverting (50) into (43) yields (44) once again but with 

k; C . (51) 

" 8 I D (/O 

In all of the above methods the time derivatives of ail of the components of 
(f> with the exception of the isotropic one, n, are neglected when compared to 
the streaming and diffusive terms. 
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Earl'^ has compared the three approaches extensively and has concluded 
that the third one is the correct approach to employ. While we hesitate to call 
any approximation scheme "correct," we do concur that the third method has 
many arguments in its favor and is the result of a fully systematic approximation 
scheme. 

However, it is known that if « (1 representing isotropic scatter- 

ing, all three methods lead to the identical value for ^ .| . Therefore if the form 
of is not too different from that of isotropic scattering one would not ex- 

pect it to matter very much which approach was used and convenience should 
dictate, especially in the field of astrophysics where the quantitative input is 
usually of a somewhat approximate nature. 

We have computed values of x,| by both the Legendre expansion method (k[') 
and the third or inversion method ) for the four sets of parameters for which 
we have complete curves of The results are shovT. in Table 1. 

As would be expected, the differences are most pronounced where - * 0) 
is the smallest. The case t = 1, ■,, = .1 for the perpendicular slab may also be 
compared with the results of Goldstein^^ who calculated by an independent 
non-linear method'^ and used Eqs. (45) and (51) to calculate Kjj- and (Goldstein 
uses the symbol rather than our .) Goldstein obtains » 115 and »• j| ■= 165. 
The near identity of with our results reflects the insensitivity of (cf. Eq, 45) 
to D ( = 0), while his larger « 1 is the result of the fact that his non-linear method 
leads t i smaller D than ours in the region " 0 (cf. Eq. 51). 

MM 

Sir ' *<11 depends inversely on the magnitude of a knowledge of how 
scales with the parameters and ■, would be of use in estimating the way «,| 
varies with the power in the fluctuating fields and particle rigidity. Since 
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K|| depends on the shape of as a function of /x as well as the overall mag- 
nitude a simple scaling will not tell the whole story. However, since the value of 
in Eq. (51) depends rather sensitively on the value of at /x* 0 we have 
chosen this value as the scaling quantity of interest. We have therefore computed 
DPAfix* 01 for the nine different cases resulting from the combination of the 3 e- 
values .25, 1, and 4 with the 3T?-values .05, .1, and .3, and we have attempted from 
these 9 cases to fit (m = 0) to the form 

DPA 0) - rr 

For each value of € we determined that m w'hich gave best fit to the theoretically 
determined values for the 3 v's. The linearity of log [D^(/x= 0)] vs log ry is ex- 
cellent for each value of e. The optimum m- values are 2.24, 2.30, and 2.34 for 
the cases -x * .25, 1 and 4 respectively. 

Similarly we determined q's holding each of the 3 values of x, fixed. The 
linearity of log = 0) vs log f at constant xy is again excellent. The opti- 
mum q-values a.e .58 -05), .60(»/ * .1) and .61 (•/ * .3). Combining fits we 

propose the empirical scaling 

(y; r: 0 ) (52) 

A similar effort to scale the values of obtained from the Monte Carlo 
simulations described in the companion paper yields^'^ 

Dj^; (m X 0)- (53) 

It is difficult to say whether the difference between the indices in Eqs. (52) and 
(53) is significant. But it is interesting to note that if one adopts the heuristic, 
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mirroring picture described in Section HI one comes up with a scaling 

- 0)- ^2.5 g.S 

intermediate between the two. 

V. Range of Validity of Partially Averaged Field Theory 

All theories, whether quasilinear or non-linear, compute the effects of the 
random force field on a collection of particles by summing its instantaneous 
effects along some particular particle orbit or collection of orbits until the ran- 
dom force on the particle becomes self-incoherent. In the quasilinear theory 
the orbit employed is the one the particle would have if the random part of the 
force were absent and we shall investigate how rapidly the actual orbit deviates 
in a significant way from this approximate orbit. In the following we shall con- 
tinue to limit ourselves to the slab model of Section III. This limitation is for 
the sake of illustration only and does not limit the validity of our conclusions. 

The quasilinear orbit is characterized by z = p vt, /x = const, but in fact, n 
will be changed by the random field by an amount 

( 54 ) 

where in this section b.. refers to the RMS value. If du; ^ ( O then will be- 
come of order unity in a time which we shall call the deflection time where 

Td - (0 for X . > ( O (55) 

If «<o then - will oscillate with a frequency ( - > and an amplitude of order 
ou>/ < - ) *= X,. However, due to the random scattering of the pitch angle, the deflec- 
tion in will grow in time as 
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where is the pitch angle diffusion coefficient calculated by the theory in 
question. In this case the deflection time is given by 

T - (D ■)“* for bv « ( - > (57) 

Clearly the random force must become incoherent in a time short compared 
to but there is another time of interest. The deflection time is a measure of 
how rapidly the particle orbit becomes unreliable in velocity space but the orbit 
can also become unreliable in coordinate space (in this model the z coordinate 
is the only one that matters). If is uncertain by an amount : (t) the z coordi- 

nate will be uncertain by an amount 


• z - V 
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(V) dt' 


(58) 
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If this uncertainty in coordinate space grow s to be of the order of the correlation 
length 7.^. which is the characteristic scale of variation of the random field, the 
unperturbed orbit becomes quite useless. The time in w'hich this occurs is called 
the deviation time and is given by 
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We further require that this deviation time be long compared to the time required 
for the random force to become self-incoherent. 

Just as we have defined two time scales over which the approximate particle 
orbit becomes unreliable we may define two coherence time scales, one for con- 
figuration space and one for velocity space. This is possible because the random 
force ■ F = (q/c) v x 5B (z) depends on both z and v. Clearly the shorter of these 
two times is the time over which F becomes incoherent. 

The correlation length z^ is defined as the length scale over w'hich the ran- 
dom field > B becomes self-incoherent, it is therefore natural to define the cor- 
relation time as the time required for the particle to travel one correlation 
length along its unperturbed orbit. In quasi -linear theory this is given by 

A. (GO) 

pv 

We may further define the coherence time as the time required for the velocity 
of the particle to become self-incoherent as it travels along its unperturbed orbit. 
In quasi-linear theory the orbit is fixed for all realizations of the ensemble so the 
particle's velocity has no stochastic component and hence never becomes self- 
incoherent. Therefore obviously in quasi-linear theory t-' = > and is the de- 
termining coherence scale. This is not the case in PAF theory. 

We have defined two time scales over which the orbit used in calculations 
becomes unreliable, one for configuration space and one for velocity space. In 
the same vein we have defined two time scales over which the random force lie- 
comes self-incoherent, again one for configuration space and one for velocity 
space. In Table 2 we list these four time scales for future reference. We require 
of any theory that 
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shorter of T . T ' < shr>rler of T,. , T . . 
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In Table 2 we have listed the symbol, name and meaning of each of the four time 
scales that we have defined here. 

We turn now to a comparison of the ratios / t ^, /■ ^ etc. as given by 

quasi-lineai dieory and PAF theory. In the interest of saving space we shall not 
write down all of the time scale ratios that can be formed although it is a simple 
matter to do so. We shall rather only explore the consequences of the constraint 
(61). 

In quasi-linear theory the requirement leads to the constraint 
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For strong random fields this condition cannot be met for all (if any) values of /.< 
but for weak random fields the condition can be met if the power spectrum goes 
i k“^ or steeper for high frequencies (i.e. ^ ) goes to zero with linearly 

or faster). 

If, however, one turns to the requirement that ^ ^ ^ one finds for ' 

1/2 
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not much different than before but for . < ■ > the requirement is 

> -■ < > 
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a condition that cannot be met for all values of ^ . We see therefore that quasi- 
linear theory will always have trouble as /j * 0. 

The various non-linear theories and PAF theory in particular were intro- 
duced to obviate this problem in the vicinity of ■ 0. This is achieved in PAF 
theory by allowing the fluctuating part of the field to influence the orbit to lowest 
order. If - « <' > this influence may be approximated 

z(t) - ^ <' > VI (65) 


and hence 


/ 

C 


( 66 ) 


Since this is equivalent to replacing /- with an effective , , ^ . < ■ ) 

a glance at (62) and (64) indicates that the troubles at (. •= 0 are avoided. 

However, for the case of strong fluctuati.ug fields, < > this substitution 

of , for /. is not a valid approximation. In fact for certain configurations of 
the fluctuating fields (other than the slab model) the particles can be trapped in 
the vicinity of the origin for an indefinitely loug time and • _ can become essen- 
tially infinite. In this situation we are saved by the fact that the coherence time 
r' which was infinite in quasilinear theory is finite in PAF theory and can lx.‘ the 
shorter of the correlation/coherence times. 

In PAF theory the velocity vector is turned not at the fixed frequency ; 

I'ut at the frequency characteristic of the partially averaged field 


^ I <B) ^ (67) 
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During final averaging this frequency will Ije spread over an amount /.^ where 

(■, f or < ' > 


A . ^ ) for-;. « < . > 


( 68 ) 


The velocity vector will thus become self-incoherent in a time t ^' given by 

/ ' - ( A. . ) " ' - I' . “ * for ■ - < - ) 
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We now find for the various ratios 
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Since as v - 0. a 



where U ■ ■ , the only ratio that does not go to zero as v - 0 is ' / ’,j in the strong 
field case. So we see that even in the case of strong random fields and low rigidity 
particles there exists an effect that destroys the coherence of the random force at 
least on the same time scale that the orbit becomes unreliable. 

Non-linear theories represent an improvement over quasi -linear theory lor 
one primary reason: they employ a more realistic description of particle orbits 
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in determining the effects of the fluctuating fields on the particles themselves. 

We have seen that one of the effects of this description is the limiting of the cor- 
relation time to a more realistic value such that the resulting calculation remains 
of marginal validity for a much wider range of parameters. This is true for PAF 
theory and would presumably be true for any non-linear theory that correctly de- 
scribes the loss of coherence over'the ensemble of the particle's velocity vector. 
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Conditional Averaging for Multi-Dimensional 
Gaussian Processes 


Consider a six dimensional vector x made up of two three dimensional 
vectors x and x, 

_X ■ (’la* 

We wish to find the average value of Xj^ under the condition that x^=x,y. If 
the joint probability density for x^ and x,^ is given as a function of k as P(x) we 
have 
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Note that 


-b 






0 


1 i k * X 


where k like x is six dimensional. Further, writing the delta function as a 
Fourier integral 
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>77)3 J 


k,' - 0 
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where k' is also six dimensional, we may insert these expressions into (A2) to 
obtain 




(■ 


'•--a 


■’iSi, 


d 6 X e ' ' " P ( X ) 


( 277)3 


k = k' = 0 


(A3) 




where k = k + k'. 

'~s^ 

The quantity in square brackets is the characteristic function of the distribu- 
tion function P(x) 

Characteristic fn. P l 

.31 


If X is a Gaussian process the distribution P(x) is a multivariate Gaussian' 

(A4) 


P(xJ - |Mi~3 exp 1 X • M"‘. x__j 

■^Tiere M"* is the inverse of the covariance matrix M. i.e,. M.^ <^x. x^"^ an 

obviouslj' symmetric matrix. In complete analogy with a one dimensional Gaussian 
distribution it is straightforward to show that 

(A5) 


P (r ) - exp • M * « J • 


If we now break M up into the three by three matrices S 

C = <;|^Xa x,,^, x,,^ where we have assumed that the marginal distributions 


for X and x, are identical we have 


M 


S C 


s 




and K • M • K = K, ,♦ S • K + K. • S • K, + k * * «|, + k,, * ’ k 
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Nv 
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Tims 




l’(K) 


k = k. = 0 
<^0 




K • M • ic I IM/c) 

I 


k = k^ =0 

'^D 


= - ('^1, ’ S ’’ 


>5^ 


k = k^ =0 

'“\^D 


= - k ' • C' cx 




p f- — k ' • S • k . 

\ ^ ~a ^ <^d j 


(kb) 


We have now 




f . -.k' • ^ -1 k' 

/ d ' k' e ~ e 2 
1 J ~-‘ 


S • k 


- k' • C’ 


Ctt)-’ 


P(X.) 


J_ 

3x„ ( 


1 / 1 -'k* * \ ~ , 

r /d'^k'e~" ~=P(k' 

2n)^ y ~a ~a 


mj 

/N^ 


r '* 1 

• (' 

; — 

Ox ~ 

_ ->_/U _ 



l’(x^,) 


= (X • S*' 


) 


c 

?w 


= X,, • (S-‘- O 
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APPENDIX B 
The Gyrotropy of <1^ 

Consider the equation for <^f')> 

.D(,)<f> <B1) 

at 

where D(/ ) represents the complete collision operator. If D( ?) were not a func- 
tion of the gyro-phase / , we could average this equation over and remove it 
from the problem but the t' dependence of D and <^f)> will average together to 
produce an unknown phase correlation. 

Let us break <(f)> into a gyrotropic part <(f^ ^ and a i dependent part <^f)>i. 

If we indicate the average over gyrophase of a quantity by placing a bai over it, 
we further specify that <(f^* i * 

Inserting in (Bl) and averaging over / gives 

-1^ D<f>„.57o; •’*2) 

Subtracting (B2) from (Bl) gives 

We now shall es* n;uv T terms in (B3) as follows 

- <f>,: D <f> ^ < > <f> 

to obtain 
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where we now assume n'^ « I and we neglect 77 2 <4^ </f ^ ^ compared to<w)><^f^j. 
Integrating (B4) with <^f)> ^ (t *= 0) = 0 gives 

<fOi ^ f <")> (t - oKOo 

Jo 




for t » 

We see from (B5) that <(f)>j is 0 ( 7 ]^) if is 0 ( 1 ); therefore the second Pjrrn 
on the right hand side of (B2) is of 0 ( 77 '*) and may be neglected. 


Case 



1, r,- ,05, islab 

456 

622 

1, T) *: .1, i slab 

115 

137 

1, r, = .3, Islab 

13.06 

13.48 

1, . 1, oblique slab 

162 

367 









Table 2 
Time Scales 


Symbol 

Name 

Meaning 

D 

Deflection time 

Time required for orbit to become 



unreliable in velocity space 

T 

(i 

Deviation time 

Time required for orbit to become 



unreliable in configuration space 

7 

Correlation time 

Time required for magnetic field 



to become self-incoherent 

7 * 

Coherence time 

Time required for velocity to be- 



come self-incoherent 




I 


\ 


i 

J 1 


Figure Captions 

Figure 1. The diffusion coefficient for the perpendicular slab model when 
t * / , m / </-^/p « 1 and '/ - = .05. 

Figure 2. vs. m for the perpendicular slab model when i = 1 and , = .1 

Figure 3. vs. /. for the perpendicular slab model when e = 1 and , * .3 

Figure 4. The p^, p^ plane showing regions where particle exits toward +z 

( ^ 0), particle exits toward -z ' 0), and where mirror- 

ing can possibly occur (indicated by "M?”). 

Figure 5. vs. /j for the oblique slab model when = 1 and ■, ■= .1 
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Figure 4. The p^, plane showing regions where particle exits toward -tz ( (a) ' 0), particle 

exits toward -z ( ( / ' " 0), and where mirroring can possibly occur (indicated by "M?"). 
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